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ABSTRACT 


The  rotationally  inelastic  diffraction  of  from  a 
corrugated  surface  is  investigated  using  a  mean  trajectory  model. 
The  center  of  mass  motion  is  treated  using  Gaussian  wave  packets, 
which  propagate  on  a  rotationally  averaged  potential.  This 
trajectory  in  turn  drives  the  rotational  transitions.  The  method 
is  non-perturbative  and  allows  for  changes  in  m ^ ,  the  rotational 
orientation,  and  agress  well  with  recent  close  coupling 
calculations.  A  connection  is  also  made  with  the  recent 
semiclassical  trajectory  work  of  DePristo.  The  effects  of  the 
attractive  well  depth  on  rotational  excitation  and  diffraction  are 
considered. 
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I .  INTRODUCTION 

The  interest  in  the  collision  of  the  lightest  diatomics 
with  surfaces  has  been  increased  by  recent  developments  in  the 
experimental  techniques.^  ®  By  using  resonant  multipohoton 
ionization,  Zare^  and  his  colleagues  have  shown  that  it  is 
possible  to  measure  the  rotational  distribution  of  molecules 
scattered  from  a  surface.  It  is  now  reasonable  to  hope  that  a 
combination  of  this  method  with  the  existing  time  of  flight 
measurements  may  permit  us  to  study  the  collision  of  light 
diatomics  with  surfaces  in  much  greater  detail  than  possible 
until  now. 

Such  developments  will  stimulate  and  help  theoretical  work 
that  intends  to  establish  a  connection  between  various  features 
of  the  molecule-surface  potential  and  the  observed  scattering 
behavior,  and  to  explain  how  the  coupling  between  the  center  of 
mass  motion  (undergoing  diffraction)  and  the  rotational  motion 
(undergoing  rotational  excitation  and  rotational  polarization) 
influences  the  outcome  of  the  scattering  events.  From  a  more 
practical  point  of  view  these  phenomena  may  become  useful  surface 
probes  since  they  are  sensitive  to  both  the  long  range  (i.e., 
through  diffraction)  and  the  short  range  (i.e.,  through 
rotational  depolarization  and  excitation)  corrugation  of  the 
surface.  They  can  also  provide  a  good  testing  ground  for  various 
approximate  theories  of  quantum  dynamics  and/or  potentials. 

In  spite  of  the  existence  of  many  interesting  theoretical 
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papers  there  is  still  a  need  for  finding  simple,  efficient, 

reliable  and  illuminating  approximate  methods.  The  earliest 
work,  ®  used  the  coupled  channel  approach®^’  which  is  exact  in 
principle,  but  is  limited  in  practice  by  the  dramatic  increase  in 
the  number  of  channels  with  the  kinetic  energy  of  the  incident 
particles.  Furthermore,  the  complexity  of  the  coupled  channel 
theory  for  scattering  by  rigid  surfaces  is  such  that  one  cannot 
be  optimistic  about  the  possibility  of  Incorporating  lattice 
motion  without  using  perturbation  theory. 

Subsequent  work,  searching  for  both  computational  simplicity 


and  an  intuitive  picture  of  collision  dynamics,  used  a  variety  of 
approximations . 

g 

Garibaldi  et.  al . ,  considered  a  hard,  corrugated  surface 
and  used  Rayleigh's  approximation  to  deal  with  diffraction.  This 
method  eliminates  many  of  the  interesting  features  of  the 
potential  (i.e.,  it's  "softness"  and  the  attractive  well)  and 
assumes  that  the  curvature  and  the  amplitude  of  the  corrugation 
are  negligible  with  respect  to  the  wave  length  associated  with  the 
center  of  mass  motion.  Furthermore  it  does  not  allow  iim.  x  0 
transitions  which  change  the  projection  of  the  angular  momentum  on 
the  normal  to  the  surface.  Thus  it  cannot  treat  correctly  the 
rotational  polarization. 

Several  groups  have  used  the  sudden  approximation,^*^  which 

assumes  that  the  transfer  of  parallel  momentum  takes  place 

suddenly  at  the  point  where  the  incident  molecule  hits  the 

surface.  As  a  result,  the  variables  X  and  Y,  which  give  the 

location  of  the  center  of  mass  in  a  plane  parallel  to  the  surface, 

can  be  treated  as  parameters  and  the  S-matrix  S.  ,(X,Y)  (j  and 

m  are  rotational  quantum  numbers)  can  be  c<.ilculated  for  each  value 

of  the  pair  (X,Y).  Diffraction  is  then  introduced  in  the  theory 

by  expanding  S.  ,(X,Y)  in  a  Fourier  series  in  which  the  wave 

jm-*j  m 

vectors  are  reciprocal  lattice  vectors;  the  S-matrix  elements  for 
various  diffraction  peaks  are  given  by  the  Fourier  coefficients  in 
this  series.  Since  the  sudden  approximation  is  an  exponentiated 
perturbation  theory,  it  is  expected  to  be  useful  for  low 
corrugation  only. 

The  semi-classical  perturbation  theory  of  Miller  and 
1 3 

Smith  was  applied  to  diatomic  scattering  from  surfaces  by 

1 4 

Hubbard  and  Miller.  They  used  a  zeroth  order  classical 
trajectory  in  which  the  periodicity  of  the  surface  is  introduced 
in  the  motion  parallel  to  the  surface  by  adjusting  the  trajectory 
to  emerge  from  the  collision  in  one  of  the  directions  required  by 
the  kinematic  restrictions  due  to  parallel  momentum  conservation. 
These  trajectories  are  used  to  calculate  perturbat ively  the  phase 
change  caused  by  corrugation  and  rotational  excitation.  This 


elegant  theory  provides  relatively  simple  analytic  formulae  for 

the  S-matrix  and  it  can  be  modified  to  take  into  account,  by 

1 5 

perturbation  theory,  the  effect  of  lattice  motion.  The 

disadvantage  is  that  one  expects  the  theory  to  work  only  for  small 

corrugation  and  when  rotational  excitation  probabilities  are 

small.  Furthermore  much  of  the  simplicity  is  gained  by  using  the 

potential  proposed  by  Lennard-Jones  and  Devonshire^®  which  does 

not  depend  on  the  azimuthal  angle.  Such  a  potential  does  not  lead 

to  correct  rotational  polarization. 

1 7 

Finally,  Whaley  and  Light  have  used  a  distorted  wave 
approximation.  The  zeroth  order  wave  function  is  obtained  by  the 
coupled  channel  method  and  describes  the  rotational  excitations 
caused  by  the  zero  wave  vector  Fourier  component  of  the  potential 
(i.e.,  this  is  a  "flat"  potential  obtained  by  averaging  the  true 
potential  in  the  directions  parallel  to  the  surface).  The  other 
Fourier  components  are  used  to  generate,  by  perturbation  theory, 
diffraction  and  rotational  polarization  changes  (i.e.,  Am?iO 
transitions ) . 

In  the  present  paper  we  develop  a  non-perturbative  method 

based  on  using  Gaussian  wave  packets  (GWP)  to  treat  the  center  of 

mass  motion  and  the  coupled  channel  method  to  deal  with  rotations. 

If  the  rotational  excitations  are  forbidden  throughout  the 

1  8 

collision,  the  theory  reduces  to  the  Drolshagen-Hel ler  theory  of 

atom  surface  diffraction.  If  the  coherence  between  the  scattered 

packets  is  suppressed  the  theory  reduces  to  an  improved  version  of 
1 9 

DePristo's  theory  of  rotational  excitation. 

The  method  presented  here  is  the  simplest  among  several 
versions  developed  by  us.  It  keeps  the  lattice  rigid  and  it  uses 
what  we  call  a  mean  trajectory  approximation.  The  latter  treats 
in  an  average  way  the  manner  in  which  the  energy  exchange  between 
translation  and  rotation  affects  the  translational  motion;  for 
this  reason  effects  such  as  surface  trapping  induced  by  rotational 
excitation,,  or  selective  adsorption  are  not  contained  in  the 
theory:  furthermore  the  post-collision  translational  energy 
distribution  is  not  given  accurately.  An  improved  method  which 


associates  one  GWP  with  each  rotational  excitation  and  in  which 

the  motion  of  the  lattice  is  incorporated  by  using  a  Langevin 

equation  is  now  being  implemented. 

Our  hope  is  that  if  successfully  completed  this  program 

will  provide  an  efficient  method  which  treats  non-perturbatively 

the  effects  of  lattice  motion  and  corrugation  to  yield 

diffraction  and  rotational  excitation  and  polarization,  and  can 

also  deal  with  situations  in  which  the  surface  is  disordered. 

The  main  criticism  of  this  method  is  that  the  use  of  GWP ' s  to 

solve  time  dependent  quantum  problems  is  not  a  well  defined 

approximation  in  which  estimates  for  the  error  and  the  range  of 

validity  of  the  theory  are  provided.  Thus  confidence  in  the 

reliability  of  the  method  can  be  established  only  by  extensive 

20 

numerical  studies.  However,  the  results  of  the  present 

calculation  for  H.  diffraction  and  rotational  excitation  compare 

^  8b 

well  with  the  exact  close  coupled  results,  as  well  as  with  the 

*1  17 

other  existing  approximate  methods*  '  in  the  parameter  range 
where  the  latter  are  reliable. 


II.  THEORY 

11. 1  Introductory  Remarks 

The  idea  of  using  Gaussian  wave  functions  to  calculate  time 

dependent  quantum  processes  was  pioneered  by  Keller  and  his 
2  1 

coworkers.  The  method  has  great  flexibility  in  regards  to  the 

manner  in  which  the  wave  function  is  constructed  and  in  which  it 

is  propagated.  Due  to  the  complexity  of  the  present  problem  we 

are  using  in  this  paper  the  simplest  choices  available.  In 

assuming  the  form  of  the  wave  function  we  use  a  hybrid  which 

treats  the  rotational  motion  by  the  coupled  channel  method  and  the 

translation  of  the  center  of  mass  by  using  a  GWP .  In  combining 

2  2 

the  two  we  use  the  mean  trajectory  approximation  (MTA)  which 
takes  into  account  the  effect  of  the  rotational  excitation  on  the 
motion  of  the  center  of  mass  in  an  average  way.  We  expect  MTA  to 
be  satisfactory  for  experiments  using  poor  angular  and  time  of 
flight  resolution.  An  improved  wave  function  form  is  provided  by 
the  multiple  trajectory  wave  function  which  associates  a  GWP 
(describing  the  translational  motion)  with  each  rotational 
transition. 

The  propagation  method  used  here  is  the  one  we  called 
20 

previously  the  simplest  Heller  method  (SHM)  which  assumes  that 

the  GWP ' s  are  decoupled  and  stay  narrow  througnout  the  collision. 

Even  though  these  conditions  are  not  well  fulfilled  SHM  works  well 

for  atom  diffraction.  Therefore  the  GWP  '  s  are  used  here  for 

the  diffraction  part  only;  SHM  should  be  satisfactory  in  the 

present  case  as  well.  The  propagation  scheme  can  be  improved  by 

using  the  minimum  error  method  in  which  the  Gaussians  are  coupled 

20 

and  can  have  arbitrary  widths. 

11 . 2  The  Hamiltonian 

We  use  a  coordinate  system  with  the  X,Y  axes  fixed  and 
contained  in  the  surface  plane  and  the  Z  axis  pointing  towards  the 
vacuum.  The  nuclei  of  the  diatomic  are  located  at  and  and 
the  center  of  mass  is  at  The  interatomic  distance  is  held 

fixed  (the  rigid  rotor  approximation)  since  the  translational 
rotational  and  phonon  energies  are  too  small  to  cause  v :  bra  t  i  ot.a  , 


excitations.  The  rotational  motion  is  described  by  the  polar 
angle  ©  (between  r  and  OZ)  and  the  azymuthai  angle  o  (between 
projection  of  r  in  the  XOY  plane  and  the  OX  axis) .  The 
Hamiltonian  is 


I 


H  =  -{h^/2m)  vj  +  (2u|rr)  t  +  V{r^.r^) 
R 


( II  .  1) 


with  m  =  m^+m^,  r  =  u=m^m2/ (m  j^+m^ )  •  ^  is  the  angu..ar 

momentum  operator.  The  potential  is 


V{r^,e^)  =  Z  D . {exp[-2a. (z.-z^) ]-  2©expr -a . (z . -z^ ) ] ) 

(II. 2) 

-p^D.  exp[-2a^(z^-z^)  ]  [  cos  (  2ttx^/c ^  ) +cos  (  2TTy  .  / c^  )  ]  • 

Here  r.  5  (x.,  y.,  z.)  and  c  and  c„  are  the  lattice  constants. 
When  the  parameter  e  is  zero  the  potential  has  no  attractive  well 
9=1  introduces  a  finite  well  depth. 
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11,3  The  mean  trajectory  wave  function 

The  mean  trajectory  approximation  (MTA)'"''  uses  a  wave 
function  of  the  form 


N  n 

iii(^,e,<ii;t)=  Z  CG{^;t)  Z  c.(t)  y,  (  9  ,  i  )  exp{ -i  .  t  h 
,  a  a  .  ,  oti  j.  1 

a=l  1=1 

(  II  .  3 


The  index  i  is  a  symbol  for  (j,m)  which  labels  the  spherical 

harmonics  Y.  ( © , p )  s  Y . { e , d ) .  These  satisfy  the  equation 
j  m  1 


(2u|r!^)"^  L  ^  Y.  (©,4>)  =  (j(j  +  l)h^/2u|r|^)Y  (©,o;=e  Y  ie,0) 

J  m  J  m  J  J 


•1 


(11.41 


The  Gaussian  wave  packet  G  is  defined  bv 

a 


G  (R;t)  =  exp{ (i/h) [ (R-R  (t))-  A  ( t ) • ( R-R  ( t ) ) +P  ( t ) • ( R-R  (t) 
a  a  0£  ot  ot  ot 


II  .  5  ’ 


The  matrix  A^(t),  called  in  what  follows  the  width  matrix  of  the 

packet  a,  has  complex  elements;  it's  imaginary  part  controls  the 

extent  of  the  packet's  spatial  localization;  the  real  part 

contributes  to  the  phase.  The  usual  definition  of  a  GWP  includes 

a  complex  phase  fator  y^(t).  We  prefer  to  incorporate  this  factor 

into  the  coefficients  c  (t).  The  constants  C  are  chosen  so  that 

a  ^  a. 

at  time  t=0  the  sum  EGG  (R;0)  is  a  plane  wave  in  the  coordinates 

a  a 

X  and  Y  and  a  normalized  GWP  in  the  direction  OZ .  The  choice  of 
the  initial  wave  function  is  discussed  in  detail  in  Sections  II. 5 
and  I I I . 1 . A . 

The  present  combination  of  a  GWP  representation  for 
translation  and  a  coupled  channel  representation  for  rotation  is  a 
new  application  of  the  Gaussian  wave  packet  method,  which  we  have 

found  useful  in  dealing  with  other  problems  such  as  curve 

2  33  2  5 

crossings  or  vibrational  excitation  by  collisions.  There  are 

two  reasons  for  using  such  a  mixed  approach:  the  translational 

energy  and  that  of  the  phonons  are  of  the  same  order  with  the 

rotational  spacing  and  therefore  we  need  only  few  spherical 

harmonics  to  describe  rotational  excitation;  since  the 

representation  of  the  rotational  motion  by  GWP ' s  require  a  large 

number  of  packets  it  is  more  practical  to  use  spherical  harmonics. 

On  the  other  hand,  the  use  of  GWP ' s  to  represent  the  center  of 

mass  motion  in  a  diffraction  problem  is  more  economical  than  the 

18  2  4 

customary  Fourier  series  expansion. 

II.  4  The  equations  of  motion  and  their  interpretation 


1 1 . 4 . a .  The  equations  of  motion 

To  obtain  the  time  evolution  of  G  (^;t)  and  c  (t)  we 

ot  Of 

introduce  the  wave  function  (II. 3)  in  the  time  dependent 
Schrodinger  equation  and  operate  from  the  left  with 


Jyy-Ti 


(II. 6) 


2tt  tt  , 

/  dd>  J"  de  sine  Y^^(e,<J)).  This  leads  to 


a  R  tt 

=  Z  G^(^,t)  Z  c^.(t) 

tt  i 

where 


2TT  TT  * 

Vj^i(^)=  Xd<b  ;  sineds  Yj^(e,4>)  V{R,e,4>)  Y^(e,4>).  (II. 7) 

r\  o 


To  solve  equation  (II. 6)  we  use  a  method  proposed  by  Heller 

which  replaces  the  potential  appearing  in  Eq.  (II.  6)  with 

a  second  order  Taylor  exoansion  in  cowers  of  R-R  (t).  To  obtain 

“  ce 

differential  equations  for  c^^{t)  and  the  parameters  appearing  in 
the  GWP ' s  we  set  to  zero  the  coefficients  of  various  powers  of  R- 
^^(t).  These  manipulations  lead  to 


=  -(2/m)  (t)  •  A^  (t)  -  (1/2)  K  ,  (II. 8a) 


?^(t)  =  (1/m)?  (t 

tt  tt 


11.8b) 


?^(t)  =  c<^l^(t)c^.  (t)exp{-i(e.-ej^)t/h}3V.j^{?^)/3R^  (II. Sc) 

and 

^  (II. 8d) 


B  (t)  =  (2ra)  P  (t)-P  (t)  -  P  (t)-R  (t)  -  ihTr  A  (t)/m, 
oc  a  o:  ot  a  oc 


(11. 9a) 


and  the  components  of  the  matrix  K  in  Eq.  (11.8a)  are  given  by 


Kj^y(t)  =22  Cj^(t)c.(t)  exp{-i(e^-ej^)t/h)  (3^V_(R^)/3X^3y^) 


k  i 


ij 


II . 9b) 


with  similar  expressions  for  K  ,  etc. 


The  equations  (II. 8)  permit  us  to  calculate  the  properties 
of  a  wave  packet,  representing  the  center  of  mass  motion,  coupled 
to  the  rotational  excitations  of  the  diatomic.  The  states 
represented  by  one  term  in  Eq  (II. 3)  (i.e.  G^(R;t)  Z 
c^^  ( t )  (  © ,  4> )  exp { -ie  ^t /h)  )  have  physical  reality  only  if  a 
special  kind  of  experiment  is  performed.  This  assumes,  in  a 
highly  idealized  fashion,  that  we  have  at  our  disposal  a  beam  of 
probing  particles  which  can  interact  with  the  incident  diatomics 

3 

in  a  volume  of  order  lA  . 

Whenever  a  probe  particle  is  deflected  by  a  collision  with 
a  diatomic  the  plane  wave  representing  the  motion  of  the  center 
of  mass  of  the  diatomic  has  been  "reduced"  to  a  wave  packet  of 
linear  dimensions  of  ~  lA.  A  coincidence  experiment,  in  which 
the  arrival  at  the  detector  of  a  diatomic  scattered  by  the 
surface  is  recorded  only  if  this  event  is  simultaneous  with  the 
arrival  of  a  deflected  probe  particle,  will  measure  the 
scattering  of  the  wave  packet  prepared  by  the  collision  between 
the  probe  particle  and  the  diatomic.  Such  a  measurement  will  not 
detect  any  diffraction  since  the  packet  acts  as  a  "corpuscule" 
which  has  lost  its  wave-like  properties. 

In  the  present  diffraction  calculation  the  packets  are  used 
as  a  labor  saving  device;  we  write  the  incident  plane  wave  as  a 
sum  of  packets  states  G^,  we  propagate  these .  packets 
independently  by  using  the  Eqs  (II.  8)  and  reconstruct  the  post¬ 
collision  state  by  adding  up  coherently  the  scattered  packets. 

In  a  diffraction  experiment  carried  out  in  the  absence  of  the 
probing  beam  the  packets  have  no  reality,  in  the  sense  that  we 
are  set  up  to  detect  only  the  state  represented  by  their  coherent 
sum 

1 1 . 4  .  c  .  The  effect  of  rotational  excitation  on  diffraction. 

To  understand  how  diffraction  is  affected  by  the  rotational 
degree  of  freedom  in  the  present  theory,  it  is  useful  to  begin 
with  the  limiting  case  in  which  the  initial  state  i=l  has  j=0  and 
m=0  and  the  diatomic  remains  in  it  throughout  the  collision. 


0  for  i >  1  . 


Thus,  c.  (t)  =  exD{iv  (t)/h),  and  c.  (t)  = 

la  a  la 

Introducing  these  values  of  c .  (t)  in  the  Eos.  (II.  3)  leads  to 
the  following  simplified  equations  of  motion;  the  time  evolutio 
of  the  phase  is  given  by 


=  B^(ti  (I^(t)) 

with 

2Tr  IT 

Vii  (Rjjj(t))  =  J*  d<{)  ;  sinede  V  ( R .  e ,  <|> )  (  4tt) 
the  equation  of  motion  fSr  momentum  becomes 


(II. 10) 


(  II  .  11  ) 


P  (t)  =  -3V  (R  ),/3^  .  (11.12) 

a  1 1  a  a 

and  the  matrix  K  has  the  elements 


These  equations  are  identical  with  the  GWP  equations  for  atom 
diffraction  in  which  the  angle  averaged  diatomic  -  surface 
potential  plays  the  role  of  an  atomic-like  potential  for  the 
center  of  mass. 

A  slightly  more  complex  limiting  case  is  that  in  which  we 
consider  two  states  only  (e.g.  i=l  is  j=0,  m=0  and  i=2  is  j=l, 
m=0 ) .  It  is  of  course  a  very  poor  approximation  to  neglect  the 
states  j=l,  m=l  and  j=l,  m=-l,  etc.,  but  we  do  this  for  the  sake 
of  simplicity:  the  insights  gained  in  this  way  can  be  extended 
easily  to  more  than  two  states.  The  equations  (Il.Sd)  can  be 
written  in  the  form 


■2a' 


(11.14) 


wi  th 


B^{t)  + 


V,,  (^(t)) 


(11.15) 


=Vi2(^^(t))  exp{-i(e2-ffi)t/h)  •  (11.16) 

An  equation  for  can  be  obtained  by  interchanging  1  and  2  in 

Eq.  (11.14).  We  can  now  integrate  Eq.  (11.14)  and  obtain 


c  .  ( t ; 

la 


=  exp{(i/h)  /  CB^(t)  +  V.^{^^(T) ) ]dT}  a^^(t) 


11,17) 


with 


=  c,^{o)-{i/h)  X  dt'  )  )a^^(t' ) 


(11.18 


-  exp{-(i/h)  X  dT[e2+V22(\{T)  )-e^-V^^(R^(T)  )])  . 
o 

Using  these  expressions  we  can  rewrite’ the  final  wave 
function  Eq.  (11.3)  in  the  form 


a  t 

Hi(^,9,<t>;t)  =  Z  C  G  (^;t)  exp{(i/fi)  X  [B  (T)dT] 


(11.19) 


2  t 

Z  exp{(i/h)  X  dTV.  .  (^^(T)  )  )  (e,<t>)  . 

i  =  l  o 

Again  for  the  sake  of  simplicity  we  examine  first  the  case  of  a 
rotationally  selected  diffraction  measurement  in  which  we  record 
only  the  angular  distribution  of  the  diatomics  arriving  at  the 
detector  in  state  1.  The  scattered  intensity  for  this  process  is 
given  by 


/  fin 


X 


I  J*dRe 


X  Y^(6,<t))  sineded*  iji  (  R,  e ,  (J> ;  t )  |  = 


Z  Z  C 
a  ^ 


* 

a 


G^XG  I 
^  a ' 


|a 


la ' 


ip 


•exp{-(i/h)  XdT  [B*  +  V..(^^(t))  -  B^-V, , (^^(t) ) ] }  . 

a  11a  BllB 

o 

(11.20) 

Since  the  particle  starts  and  ends  the  collision  in  the  state  1 
we  are  considering  an  elastic  process;  however,  we  do  not  prevent 
as  we  did  in  the  earlier  discussion  "virtual"  rotational 
excitations  from  taking  place.  Since  diffraction  is  mainly  ■ 
determined  by  the  relationship  between  the  phases  of  various 
packets  appearing  in  the  sum  in  Eq,  (11.20)  it  is  useful 
to  note  that  there  are  two  kinds  of  phases.  The  phase 
X  dT[B  +V,  ,  (^  (t))]  deoend's  on  how  the  packet  G  moves  on  the 
rotationally  averaged  potential  This  is  the  phase 

that  appears  in  the  atom  diffraction  limit  discussed  above.  The 
phase  appears  because  the  diatomic  undergoes  "virtual" 

transitions  l-»2  and  2-*’l  as  it  evolves  from  the  Initial  state  1  to 
the  final  state  1.  This  phase  is  absent  in  the  atom-like  limit 
in  which  such  virtual  transitions  are  forbidden.  The  presence  of 
leads  to  a  dephasing  similar  to  the  Debye-Waller  effect  of 
the  phonons:  they  arise  because  virtual  transitions  take  place 
during  elastic  scattering.  We  expect  this  dephasing  to  diminish 
the  peak  intensity  but  not  to  broaden  the  peak  since  it  is 
unlikely  that  the  virtual  excitations  occurring  at  different 
places  on  the  surface  are  correlated. 

We  can  analyse  in  a  similar  manner  an  angular  distribution 
measurement  that  records  only  diatomics  reaching  the  detector  in 
state  2.  A  measurement  that  records  the  diatomic  regardless  of 
its  rotational  state  consists  of  a  coherent  superposition  of  the 
1-*1  and  l-*2  processes  and  therefore  consists  of  their  sum  and  the 
interference  between  them. 


While  this  analysis  is  not  quantitative  it  does  reveal  the 
basic  physical  processes  by  which  the  rotational  internal  degrees 
of  freedom  affect  diffraction.  A  very  similar  discussion  could 
be  made  if  we  were  to  consider  the  diffraction  of  a  diatomic 
which  can  be  vibrationally  excited  by  the  collision  (of  course 
the  kinetic  energy  would  have  to  be  rather  large  for  this  process 
to  occur)  or  the  diffraction  of  an  atom  (e.g.,  Li)  which  can  be 
ionized  during  the  collision  process. 


gl£ 


The  two  level  limiting  case  summarized  by  Eqs .( II . 14-16 ) 
can  also  be  used  to  discuss  qualitatively  the  rotational 
excitation  process.  This  process  is  driven  by  the  time 
dependence  of  the  matrix  element  V )  )  which  is  controlled 
by  the  time  dependence  of  the  center  of  mass  position  ^^(t). 

This  is  very  similar  to  the  way  in  which  the  nuclear  motion 
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induces  charge  transfer  in  two  state  models.  '  However,  the 

trajectory  ^^(t)  appearing  here  is  not  a  classical  trajectory  but 

it  contains  some  feed-back  regarding  the  effect  of  the 

instantaneous  state  of  rotational-  excitation  on  the  center  of 

mass  motion;  this  is  achieved  through  the  presence  of  the 

potential  Z  Z  *^10^^^  Eqs.  (11.8).  Since 

this  feature  is  typical  of  MTA  we  are  not  discussing  it  further 

and  refer  to  our  previous  work  regarding  the  use  of  MTA  in  curve 

22a 

crossing  problems. 

11.5  The  construction  of  the  initial  wave  function. 


The  construction  of  the  initial  wave  function  from  an 

expression  of  the  form  (II. 3)  follows  the  work  of  Drolshagen  and 
18  2  4 

Heller  and  Jackson  and  Metlu.  Since  we  intend  to  describe 
diffraction  we  must  insure  that  Oi  (^,  e ,  {J> ;  t=0 )  is  a  planar  wave  in 
the  (X,Y)  variables,  (or  that  it  represents  the  coherence 
properties  of  the  beam  used  in  experiments). 


iii(I,e,<t):t=0)  =  G(z)  Z  exp(  (  i/h)  [  (  [X-X  (0))^  + 

a=l  “ 

(Y-y^(0))^]A  +  Pj^(O)  (X-X^(O)  )  + 


with 


Py(0)  (Y-Y^(O)  )  ]  )  •  (TTli/2A) 
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G(z)  =  exp{  (  i/h)  [A(Z-Z(0)  )'‘+P  (0)  (Z-Z(O)  )  +  P  (0)2(0)]) 

z  z 


(11.22) 


The  sum  in  Eq.  (11.21)  is  the  two  dimensional  expression  obtained 
by  discretizing  the  identity 


exp{iP-R/h}  =  C  J-dR  exp{  (i/h)  t  (R-^  )•  A  -(^-R  ) 

o  o  o 

?• 

o  o 


(II .  23.) 


(C  is  a  normalization  constant)  in  order  to  build  a  two  dimension 

plane  wave.  The  positions  X^,  are  taken  on  a  grid  covering  the 

unit  cell.  components  of  the  mean 

parallel  momentum  of  the  incident  planar  wave.  The  factor 
3/2 

(Trh/2A)  normalizes  the  Gaussians.  The  Z  dependence  of  the  wave 

function  is  given  by  a  GWP ;  the  manner  in  which  this  choice 

affects  beam  coherence  has  been  discussed  by  Jackson  and 
24 

Metiu.  The  position  Z(0)  is  taken  outside  the  surface-molecule 
interaction  region.  This  procedure  fixes  the  constant  C 
appearing  in  Eq.  (II.  3)  to 

=  exp[  ( i/fi)P^.R^]  (7Th/2A)^''^  .  (11.24) 

Thus  all  the  initial  conditions  for  the  equations  of  motion  (II. 8) 
are  specified  as  long  as  numerical  values  are  chosen  for  A.  The 
choice  of  the  latter  is  discussed  in  Section  IlI.l.A. 

I I . 6  A  critical  discussion  of  the  approximation. 

The  theory  outlined  here  makes  three  approximations;  The 
mean  trajectory  approximation  imposing  a  wave  function  of  the  form 
(II. 3-4);  it  assumes  that  the  Gaussians  can  be  propagated 
independently  and  the  potential  expanded  around  the  center  (SKM); 
and  it  uses  a  rigid  lattice.  We  discuss  below  the  limitations  of 
these  assumptions,  and  outline  our  current  work  aimed  at  improving 
the  theory. 

As  we  have  discussed  in  Section  II. 2  in  the  limit  of  elastic 
scattering  (i.e.,  ~  I  =  0.  n  =  2,  ....), 

the  present  calculation  is  identical  to  that  of  an  atom  scattering 
by  a  surface.  Numerical  work  has  shown^®'^^  that  SHM  with  a  wave 


function  of  the  form  (II. 3-4)  gives  a  satisfactory  description  of 

the  problem  even  though  some  questions  remain  concerning  the 
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ability  to  describe  the  coherence  of  the  incident  state.  The 

wave  function  (II. 3-4)  leaves  out  some  important  details  on  the 

description  of  the  rotational  excitation's  effect  on  the  motion  of 

the  center  of  mass.  The  rotational  spacing  for  is  of  the  same 

order  of  magnitude  as  the  translational  kinetic  energy  and  the 

rotational  excitations  cause  the  appearance  of  new  diffraction 

peaks.  Their  positions  can  be  obtained  from  energy  and  parallel 

momentum  conservation  conditions.  The  mean  trajectory 

approximation  is  not  likely  to  give  the  correct  intensity  of  these 

peaks.  To  understand  why  we  compare  the  MTA  wave  function  (II.3- 

2  3 

4)  with  the  multiple  trajectory  approximation  form 

♦  =  z  Z  (11.25) 

Ot  i 


which  associates  one  GWP  with  each  rotational  state  rather  than 

one  GWP  for  all  of  them.  If  we  consider,  for  simplicity,  a  two 

state  problem  (i.e.  two  rotational  states  and  Y^),  the 

collision  described  by  the  wave  function  (11.25)  can  generate  two 

packets  for  each  a;  one,  G  , ,  having  the  same  energy  as  the 

incident  packet  and  the  other,  '  having  the  incident  energy 

minus  the  energy  of  the  rotational  excitation.  The  coherent  sum  I 

G^2  used  to  calculate  the  intensity  of  the  new,  rotationally 

shifted  diffraction  peaks  while  Z  G  ,  generates  the  intensity  of 

the  rotationally  elastic  diffraction  peaks.  The  behavior  of  the 

single  packet  G^  used  by  MTA  is  an  average  of  the  behavior  of 

G  ,  and  G  „ . 
al  a2 

The  usefulness  of  MTA  depends  on  the  measurements  being 
analyzed.  If  one  is  interested  in  rotational  population  measured 
with  poor  angular  resolution  MTA  should  be  satisfactory,  the 
averaging  effect  of  the  large  detector  aperture  is  simulated  by 
the  MTA  wave  function.  If  one  is  interested  in  detailed,  highly 
resolved  angular  distributions  one  should  use  the  multiple 


trajectory  approximation.  In  the  same  vein,  detailed  time  of 
flight  measurements  (TOF)  require  the  use  of  multiple 
trajectories,  while  MTA  is  satisfactory  if  the  TOF  resolution  is 
poor . 

The  SHM  assumptions  that  the  Gaussians  can  be  propagated 

independently  and  that  the  potential  can  be  replaced  by  a  local 

quadratic  expansion  around  the  center  of  the  packet,  can  both  be 

removed.  The  general  procedure  for  accomplishing  this  was 

20 

discussed  elsewhere,  and  we  hope  to  apply  it  to  the  present 
problem  soon. 

Finally,  the  thermal  motion  of  the  lattice  can  be 

incorporated  into  the  present  calculation  in  a  manner  similar  to 

.  2  0 

that  employed  in  our  work  on  diffraction  which  included  the 
classical  motion  of  the  lattice  atom  through  a  Langevin  equation. 


L 
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III.  NUMERICAL  RESULTS 

111.1.  Computational  Details 

The  numerical  calculations  reported  here  use  the  wave 

function  (II. 3)  with  the  coefficients  C  given  by  (11.24).  The 

0: 

time  evolution  of  the  wave  function  given  by  the  time  dependence 
of  the  parameters^  (t),  ^  (t),  ^  (t)  and  c  (t)  is  determined  bv 
solving  the  first  order  differential  equations  (II. 8).  The 
choice  of  the  initial  conditions  for  these  parameters  and  of  the 
molecule-surface  interactions  are  discussed  below. 

111.1.  A  The  choice  of-  the  initial  wave  function 

The  proceudre  used  for  choosing  the  initial  parameters  in 
the  wave  function  was  explained  in  Section  II.  5.  We  give  some 
additional  details.  The  surface  is  a  square  lattice  with  the 
lattice  constants  c^*4A  and  C2=2A.  To  construct  the  incident  two 
dimensional  plane  wave  (see  Section  II.  3)  we  place  in  a  plane 
parallel  to  the  surface  and  located  at  a  distance  Z(0)  a 
rectangle  of  the  size  of  the  unit  cell.  We  divide  it  into  N 

CS 

subcells.  In  the  center  of  the  Gaussian  G  we  oalce  a  subcell  a. 

a 

This  fixes  X  (0),  Y  (0)  and  Z  (0).  We  test  the  results  and  see 
a  a  ot 

that  the  results  are  not  sensitive  to  these  choices  by  increasing 
Z^(0)  and/or  by  increasing  the  number  of  cells. 

The  intial  momentum  P(0)  was  chosen  so  that  P(0)  =  hk 

o 

where  k  is  the  initial  wave  vector.  The  initial  ohase  of  the 
o 

packet  G^  is  given  by  '^^(O)-P(O)  (see  Eq .  (11.23)). 

The  choice  of  the  elements  of  the  matrix  A(0)  is  more 

^  y 

difficult  and  involves  some  arbitrariness.  We  assume  that  AO'  = 

AI  where  I  is  the  unit  matrix.  Following  Drolshagen  and 
1 8 

Heller,  the  diagonal  elements  of  A  are  chosen  to  give  a  .narrow 
packet  at  the  point  of  impact  with  the  surface.  The 
justification  for  this  procedure  is  that  the  propagation  scheme 
works  best  for  narrow  Gaussians,  therefore  we  can  use  the 
flexibility  in  chosing  the  magnitude  of  A  to  make  the  packet 
narrow  at  the  most  important  point  along  the  trajectory,  i.e  . 
the  point  of  impact.  An  inexpensive  way  of  implementing  th:s 
idea  is  to  use  for  A(t)  the  free  soace  formula  and  determine  A  3 


i:. 


20 


to  give  at  the  time  of  imoact  t  the  desired  values  for  A( t  ) , 

22  ■  ^  c 

In  our  previous  work  we  have  shown  that  the  interaction  affects 

strongly  the  magnitudes  of  ImA  (t),  ImA  (t),  and  ImA  (t). 

Also,  packets  which  correspond  to  particles  of  very  small  mass 

(e.g..  He,  H^)  broaden  very  quickly.  As  a  result,  the  GWP ' s  are 

much  broader  during  the  interaction  than  would  be  desirable. 

Thus  we  are  not  achieving  the  objective  that  was  used  to  justify 

the  procedure.  Nevertheless,  the  procedure  works  in  the  sense 

that  predictions  made  for  atom  diffraction  agree  with  the  exact 

calculations,  and  the  results  are  stable  with  respect  to  small 

varitions  in  the  value  of  A{0).  Perhaps  this  happens  because 

diffraction  is  dominated  by  the  phases  of  the  GWP ' s  .and  in  SHM 

the  phases  are  essentially  given  by  the  classical  action.  Thus 

the  SHM  results  have,  roughly,  the  same  accuracy  as  those  of  the 

semi-classical  diffraction  theory.  Since  the  latter  works 

reasonably  well  so  will  SHM,  as  long  as  the  contributions  to  the 

S-matrix  due  to  the  widths  of  the  packets  are  less  important  than 

those  due  to  their  phases.  If  this  is  the  case  the  results  are 

not  sensitive  to  small  changes  in  A  (t)  caused  by  changes  in  the 

20 

choice  of  the  initial  value  A.  However  other  applicatons  of 
the  GWP  method  indicate  that  a  proper  choice  of  the  initial 
parameters  (by  varying  the  parameters  in  the  GWP  e.xpression  for 
the  wave  function  to  get  the  best  fit  of  the  initial  condition) 
is  sometimes  very  important.  Therefore  we  regard  the  freedom  in 
the  choice  of  A  as  an  unsatisfactory  feature  of  the  present 
theory . 


Ill  .  1  .B 


The  computed  quantities 


We  propagate  the  packets  for  a  time  t  sufficient  to  allow 
them  to  emerge  from  the  interaction  region.  The  observables 
computed  here  become  independent  of  the  time  when  the  packets  are 
no  longer  interacting  with  the  surface.  The  S-matri.x  for  a 
H^molecule  emerging  in  the  final  state 


q^'^^exp{  iq- R}  Y^(e,(t)', 


(III. 


I 


S'-  ={(q  )  }  I  /dR  exp{-iq  -R)  G  ( R ; r )  C  c  .  t)  ; 

nm  nm  ,  nm  ct  a  cxi 

cx=  1 

(  III  .  2  ) 

the  probability  of  reaching  the  specified  final  state  is 

i  i  2 

p  =  1 r . 
nm  nm 

The  final  momenta  must  satisfy  parallel  momentum  and  energy 
conservation  and  are  constrained  to  take  the  values 


(q  )  =  k  +  2TTn/c, 

' nm  ox  1 


( III . 3a: 


(q  )  =  k  +  2TTm/C- 

^y  nm  oy  2 


,  r  ,  .1,2 

'^z^nm  =  <^o-^^<5x'nm^ 


, 1/2 


( III . 3b) 


( III . 3c ) 


Here  1<  =  ?  /n  and  P  is  the  initial  mean  momentum  of  the  oackets 

o  o  o 

As  usual  the  integers  (n,m)  label  the  diffraction  peaks. 

The  above  equations  refer  to  scattering  from  one  unit  cell. 
The  S-matrix  for'  the  whole  surface  can  be  obtained  by  adding  the 
S-matrices  for  each  cell  weighted  with  the  the  proper  exponential 
factors  reguied  by  the  two  dimensional  peridocity.  When  this  is 
squared  it  leads  to  6  functions  which  give  (III. 3a)  and  (III. 3b). 
III. 2  Numerical  results  for  a  ootential  without  an  attractive 


In  order  to  compare  our  results  with  other  calculations  we 

use  the  potential  given  by  Eq .  (11.2)  with  e=0  and  the  parameter 

set  1  given  in  Table  I.  We  denote  this  potential  .  By  taking 

e=0,  we  suppress  the  molecule-surface  attraction.  The  incident 

molecules  move  in  a  direction  perpendicular  to  the  surface,  have  a 

kinetic  energy  of  0.1  eV  and  are  in  the  state  j=0,  m=0 

corresponding  to  para-hydrogen.  We  denote  this  calculation  . 

The  results  of  our  calculations,  together  with  those 
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obtained  by  multichannel  distorted  wave  (MDW)  ,  the  diffractive 

10c 

sudden  rotational  close  coupling  (DSCCR)  ,  and  the  exact  close 
coupled  ( CC )  methods,®'^  are  given  in  Table  II. 


The  agreement  is  rather  good,  even  over  ten  orders  of 
magnitude,  considering  that  the  mean  trajectory  method  treats  the 
translation  in  an  average  manner.  One  should  note  that  both  MOW 
and  DSCCR  are  valid  for  small  ^  and  little  off-specular 
scattering.  Thus,  our  results  should  be  similar  to  theirs  for  low 
^  systems  such  as  metal  surfaces.  When  there  is  strong  off 
specular  scattering,  which  is  the  case  for  experimental  systems 
such  as  H^/LiF  or  D^/NaF,  the  other  methods  should  break  down. 

Our  method  should  correctly  predict  the  rotational  transitions  in 
this  limit. 

Note  that  one  disagreement  between  the  present  theory,  the 
MOW  calculation , and  the  DSCCR  results  is  for  the  specular  peak. 
This  is  not  surprising  because  in  the  MDW  approach  the  zeroth 
order  state  is  taken  to  correspond  to  3=0  which  leads  only  to 
specular  scattering  and  Amj=0.  First  order  perturbation  theory  is 
then  used  to  compute  all  off  specular  scattering.  As  is  well 
known  for  all  first  order  theories,  the  zeroth  order  state  is 
unchanged  by  the  perturbations,  and  the  total  probability  is  not 
normalized.  The  problem  then,  is  that  one  can  not  calculate  the 
effects  of  having  a  finite  corrugation  on  the  rotation-diffraction 
probabilities  of  the  specular  peak,  which  is  the  only  easily 
observable  peak  for  the  low  3  systems  where  such  approaches  are 
valid , 

Thus,  we  see  discrepancies  between  the  MDW  and  exact  CC 
results  at  larger  B's  for  the  specular  peak.  The  DSCCR 
probabilities  are  very  close  to  the  exact  values,  while  the  mean 
trajectory  results  are  reasonable,  but  not  quite  as  good.  The 
three  approximate  methods  show  roughly  the  same  level  of  accuracy 
for  the  off  specular  rotationally  elastic  probabilities,  when 
compared  with  the  exact  results.  However,  the  mean  trajectory 
method  is  much  more  accurate  than  either  the  MDW  or  DSCCR 
calculations  for  off  specular  scattering  where  Amj^^O,  for  all 
values  of  3  tested.  For  the  (0,1)  peak,  the  MDW  and  DSCCR  results 
are  off  by  a  factor  of  2  or  more,  while  the  mean  trajectory  method 
gives  values  reasonably  close  to  the  exact  CC  results.  For  the 


"more  sudden"  (1,0)  diffraction  state,  this  difference  is  not  very 
large.  For  the  off  specular  |0,0>  ->  |2,0>  transition,  the  mean 
trajectory  approach  is  not  quite  as  accurate  as  the  MOW  and  DSCCR 
for  the  (1,0)  peak,  although  all  are  fairly  close  to  the  exact 
results.  The  mean  trajectory  results  are  much  better  than  both 
MDW  and  DSCCR  for  the  transition  |0,0>  -*•  (0,1)|2,0>. 

Thus,  the  mean  trajectory  method  is  a  viable  approximate 
approach  to  the  rotationally  inelastic  diffraction  problem.  Not 
only  does  it  require  a  small  fraction  of  the  computer  time  needed 

by  CC ,  MDW,  and  DSCCR  calculations,  but  it  is  easily  extended  to 

18  26 
disordered  or  finite  temperature  systems,  or  high  energy 

systems  involving  a  great  many  diffraction  states.  A  mojre 

detailed  comparison  of  the  MDW,  DSCCR,  and  CC  methods  can  be  found 

in  reference.  8b . 
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I I I . 3  The  probability  of  rotational  excitations:  a  comparison 
with  DePristo's  results. 

19 

Recently  DePristo  carried  out  a  mean  trajectory 
calculation  of  the  rotational  excitation  caused  by  collisions  with 
a  surface.  He  used  the  wave  function 


i|){e,4>;t)  =  Z  c.{t)  Y.  (6,<|))e“^i^/^ 
i  ^  ^ 


(III. 4a] 


with  the  amplitudes  given  by 


ihc  (t)  =  Z  -H  (R(t))  c 

J  ^  J  J 


( III , 4b) 


I  (R(t))  =  /  sined©  J*  d<f>  Y .  {  ©,<|> )  V{  ©,  4>  ;^(  t )  )  Y.(©,4>). 


(III. 4c) 


The  center  of  mass  motion  is  given  by 


( t )  =  -  —— 


3R(  t; 


Z  Z  c  c  H  (R(t) 
i  j  J  J 


(III . 4d) 


The  above  equations  are  a  limiting  case  of  GWP-MTA  obtained  by 

replacing  c^G^f^;t)  with  1  in  Eq.  (II. 3)  and  making  B^(t)=0  in 

(II. 9a).  A  dramatic  effect  of  the  replacement  C  G  (^;t)  with  1  is 

a  <x 

that  the  interference  between  the  trajectories  is  lost,  therefore 
the  theory  is  unable  to  give  diffraction.  For  heavy  diatomics 
this  is  not  an  important  error  but  for  such  cases  the  rotational 
spacing  is  very  small  and  the  rotational  quantum  effects  are  also 
negligible;  therefore  a  fully  classical  calculation  should  be 
adequate.  For  those  cases  in  which  the  rotational  motion  needs  to 
be  treated  quantum  mechanically  the  diffraction  effects  are  likely 
to  be  important  and  we  must  understand  how  they  influence  the 
angular  dependence  of  the  rotational  excitation  probabilities. 


DePristo's  calculation  assumes  that  this  influence  is  negligible. 

To  answer  such  questions  we  compare  a  GWP-MTA  calculation 
with  an  "incoherent  GWP-MTA",  that  is  with  a  calculation  in  which 
the  probability  P^^  that  the  particle  emerges  in  the  state  J,m  is 
given  by 

P^"  =  J-sinede;  d*  Z  |  (6,<|))  1  ^  (III. 5) 

^  o  o  o:  ^ 

instead  of 

=  J-sinede/^dcDZ  Z  e,(l>)  ]  *  (  e,(D)  ]  (111.6)  . 

jm  ^  o  a  p 

Here  .  (e,<^)  is  the  S-matrix  for  the  transition  (i=0,  m=0)  -*■ 

-  (a) 

(j,m)  at  the  angles  e,<|),  for  the  packet  G^  .  The  major 
approximation  made  by  DePristo  is,  in  essence,  the  replacement  of 
(III. 6)  by  (III. 5) . 

The  results  obtained  with  the  full  GWP-MTA  summed  over 
diffraction  peaks  (equivalent  to  Eq.  (III. 6)  and  with  its 
"incoherent"  version  (i.e.,  Eq.  (III. 5))  are  shown  in  Table  III. 
We  find  that  the  rotational  excitation  probabilities  given  by  the 
incoherent  expression,  are  systematically  larger  than  those  given 
by  the  coherent  one.  However,  the  errors  are  not  very  large  and 
therefore  we  feel  that  the  present  calculation  justifies  the  use 
of  DePristo's  method  for  light  scatterers,  when  a  qualitative 
analysis  in  intended. 

III. 4  The  effect  of  an  attractive  well 

The  recent  calculations®'^ '  of  coupled  rotational 

excitations  and  diffraction  with  which  we  compare  our  results  are 

confined  to  purely  repulsive  potentials.  The  absence  of  an 

attractive  well  simplifies  the  theory' since  it  eliminates 

selective  adsorption  and  rotationally  mediated  selective 

adsorption.  Prom  an  experimental  point  of  view  this  elimination 

2 

is  unfortunate  since  such  resonances  are  very  common. 

While  the  GWP-MTA  theory  presented  here  is  not  designed  to 
give  bound  state  resonances  it  can  be  used  to  examine  the  effect 
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of  molecule-surface  attraction  on  diffraction  and  rotational 


excitation.  In  particular  it  can  examine  whether  there  is  any 
truth  to  the  notion  that  the  effect  of  an  attractive  well  is 
roughly  equivalent  to  a  change  in  the  normal  part  of  incident 
kinetic  energy  by  an  amount  equal  to  the  well  depth.  We  find  that 
this  is  not  the  case. 

Our  first  calculation  with  an  attractive  well  uses  Eq . 

(II.  2)  with  the  parameter  set  2|p=.05  and  e=l .  We  denote  this 

potential  and  the  calculation  of  scattering  by  it  of  a  normal 

beam  with  the  incident  kinetic  energy  e^=0.1  eV  is  denoted  C^.  We 

have  chosen  Z  =1.8841A  and  D=.02  eV  so  that  we  obtain  a  well  deoth 
o 

of  0.02  eV,  and  the  potential  reduces  to  the  potential  (used 
•in  Section  III.  2)  when  we  take  e=0 .  The  results  of  the 
calculation  are  shown  in  Table  IV.  When  we  compare  them  to 
those  of  (see  Table  II)  we  find  substantial  differences.  Going 
from  to  lowers  the  intensity  of  the  rotationally  elastic 
specular  peak  and  increases  the  intensity  of  the  diffracted  beams. 
The  probability  of  the  rotational  excitation  |0,0>  -*•  |2,0>  in  the 
specular  peak  goes  down  when  the  well  is  turned  on. 

We  can  ask  whether  the  scattering  of  a  0.1  eV  beam  from  the 
potential  V^,  having  a  0.02  eV  well  is  equivalent  to  scattering 
from  the  potential  with  an  incident  kinetic  energy  of  0.1  + 

0.02  eV.  To  answer  the  question  we  compare  the  probabilities 
describing  the  latter  process  (Table  V)  to  those  of  the  former, 
given  in  Table  IV.  Clearly  the  results  are  not  similar.  Nor  are 
the  trends  in  going  from  the  calculation  to  similar  to  those 
seen  when  going  from  to  a  kinetic  energy  of  0.12  eV  on  the 
potential  V^ . 

Since  the  change  from  V^  to  V^  not  only  adds  an  attractive 
well  but  also  changes  the  slope  of  the  repulsive  wall  and  its 
effective  corrugation,  we  have  designed  a  new  potential  V  with 
the  parameter  set  3  (Table  I),  3=. 022,  and  e=l .  V^  has  a  well  of 
depth  0.02  eV  and  the  same  corrugation  amplitude  and  repulsive 
force  (for  both  e  =0.1  eV  and  e  =0.12  eV)as  the  potential  V, , 


with  p=.05.  The  calculation  denoted  C  gives  probabilities  (Table 
VI)  for  scattering  with  e  =  0.1  eV  and  from  the  potential  . 

They  are  still  rather  different  from  those  obtained  for  scattering 
by  the  purely  repulsive  potential  at  the  kinetic  energy  0.12 
eV. 

It  seems  that  for  the  situations  presented  here  the  presence 
of  an  attractive  well  has  a  significant  influence  on  scattering 
probabilities,  which  cannot  be  predicted  by  assuming  that  the  well 
merely  changes  the  incident  kinetic  energy. 
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Table  I . 

The  values  of  some  of 

the  parameters 

appearing  in  the 

1 

potential  { Eq .  II. 2). 

The  remaining 

parameters  are 

F 

1 

specified  in  the  text 

i 

Parameter  1 

2 

3 

i 

Set 

-1 

1 

a  (A  "■) 

1 . 83315 

1 . 83315 

1  .  17 

1 

2^  (A) 

0.0  A 

1 . 8841 

2 . 5 

o 

D  (eV) 

20.0 

0.02 

0.02 

1 

Cj  (A) 

4.0 

4.0 

4  ;0 

C2  (A) 

2.0 

2.0 

2 . 0 

I  J 


Table 

I 

G 

(0,0) 

(1,0) 

(0,1) 

(0,0) 

(1,0) 

(0,1) 


II.  Transition  probabilities  |j,m>  =  |0,0>  G'|j',m'> 

for  H-  on  a  rectangular  lattice  at  normal  incidence 
Results  are  for  the  gultichannel  Distorted  Wave^ 
(MDW) ,  close  coupled”  (cc),  mean  trajectory  (this 
work) ,  and  diffraction  sudden  rotation  close 
coupled"^  theories  respectively.  The  calculation 
was  performed  for  parameter  set  1  of  Table  I  for  3 
values  of  the  corrugation  parameter  p. 


t  j'm’> 


p 

|0,0> 

12. 2> 

1  2, 1> 

1  2,0> 

CC 

.005 

.904 

.832(-10) 

0 

.957(-1) 

.05 

.aeo 

.ai4(-6) 

0 

. 920( -1 ) 

.  1 

.737 

.122 (-4) 

0 

.813(-l) 

.005 

. 140(-3) 

.43l(-7) 

. 123{ -6) 

.  lll(-4) 

.05 

. 136(-1 ) 

.464(-5) 

.121{-4) 

. 109(-2) 

.  1 

.500(-l) 

.23l(-4) 

.456(--4) 

.  406( -2) 

.005 

.826(-4) 

.  .335(-6) 

. 194(-6) 

.640(-5) 

.05 

.800(-2) 

.327(-4) 

. 189(-4) 

. 624( -3) 

.  1 

.290(-l) 

.  12l(-3) 

!704(-4) 

.230(-2) 

MEAN  TRAJECTORY 

.005 

.917 

.854(-I0) 

.825(-l) 

.05 

.877 

.997(-6) 

'•.3(-16) 

.792(-l) 

.  1 

,759 

,  155 (-4) 

-.2(-16) 

.697(-l) 

.005 

. 146(-3) 

.412(-7> 

. 115{-6) 

,984( -5) 

.05 

. 136(-1 ) 

.446(-5) 

. 112(-4) 

. 930( -3) 

.  1 

,501(-1) 

. 225(-4) 

.416(-4) 

.349(-2) 

.005 

.744(-4) 

.370{-6) 

. 243 ( -6) 

.522(-5) 

.05 

.675(-2) 

.362(-4) 

. 239{ -4) 

. 526( -3) 

.  1 

. 246{-l ) 

.  135(-3) 

.911(-4) 

.  195(-2) 

Table  II.  continued 


G  ' 

p 

|0,0> 

|2,2> 

MOW 

|2,  1> 

1  2 , 0> 

.005 

.893 

0 

0 

.  106 

(0,0) 

.05 

.893 

0 

0 

.  106 

.  1 

.  893 

0 

0 

.  106 

.005 

. 138(-3) 

.467(-7) 

.  138(-6) 

. 125(-4) 

(1,0) 

.05 

. 138(-1) 

.467(-5) 

. 138(-4) 

. 125(-2) 

.  1 

.553(-l  ) 

. 187(-4) 

.552(-4) 

.500(-2) 

.005 

.898(-4) 

.634(-6) 

.465(-6) 

.932(-5) 

(0,1) 

.05 

.898(-2) 

.634(-4) 

.465(-4) 

.932{-3) 

.  1 

.359(-l) 

.254(-3) 

DSCCR 

. 186(-3) 

.373(-2) 

.005 

.905 

.843(-10) 

0 

.948(-l ) 

(0,0) 

.05 

.859 

.822(-6) 

0 

.905(-l) 

.  1 

.441 

. 103(-4) 

0 

.  157{-2) 

.005 

. 142( -3)  ■ 

.506(-7) 

.  l54(-6) 

.  117(-4) 

(1,0) 

.05 

.137(-1) 

.549(-5) 

.  149(-4) 

.  113(-2) 

.  1 

.130 

.156(-4) 

.433{-5) 

.329(-3) 

.005 

.925(-4) 

.696(-6) 

.519(-6) 

.867(-5) 

(0,1) 

.05 

.890(-2) 

. 672( -4) 

.503(-4) 

.838(-3) 

.  1 

.814(-1) 

.507 (-4) 

.  147(-4) 

.233(-3) 

Table  IV. 


Transition  probabilities  |0,0>  -►  (  n  ,  m )  |  j  '  m  '  >  for  an 
incident  beam  of  .1  eV  H.  molecules,  us ing ‘paramete 
set  2  with  ©  /  and  8=.05.  Also  included  are  the 
incoherent  ( P .:  )  and  diffraction  summed  results  iZ) 


|0,0> 

12, 0> 

|2,1> 

1  2  ,  2> 

(0,0) 

.  730 

. 237(-l ) 

^10-'® 

.282(-4) 

(±1,0) 

. 676( -1 ) 

. 197(-2) 

. 446 ( -5 ) 

.448(-4) 

(0,±1) 

.400(-l) 

. 150(-2) 

. 127(-4) 

.216{-3) 

(±1  ,±1 ) 

.333 ( -2 ) 

.  110(-3) 

. 147(-5) 

.  125 ( -4 ) 

(±2,0) 

. 121 (-2) 

.288(-4) 

. 444 ( -6 ) 

,  225 ( -5 ) 

10, ±2) 

.200(-3) 

.674(-5) 

.872(-6) 

.320{-5) 

^jm 

.951 

.  0448 

, 108(-3) 

.  152 ( -2 ) 

Z 

.  967 

.312(-1) 

,428(-4) 

.611(-3) 

Table  V. 

Transition  probabilities  |0 

,  0>,  -►  (  n ,  m)  1  j  '  , 

,  m ' > ,  and 

incoherent  and 

diffractive 

summed  results 

for  a.n 

incident  beam 

of  . 12  eV  H- 

molecules,  using  parameter 

set  1,  with  e= 

0  and  8=.057 

A 

0 

0 

1  20> 

1  21> 

1  22> 

(0,0) 

.843 

.  102 

~  1 0 

.171(-5) 

(±1.0) 

. 158(-1 ) 

. 158(-2) 

. 147(-4) 

.710(-5) 

(0,±1) 

.830(-2) 

.987{-3) 

.331 (-4) 

.591(-4) 

(±1,±1) 

.104(-3) 

.101(-4) 

<io"® 

<io‘® 

3m 

.  850 

.  148 

. 198(-3) 

, 330( -3 ) 

Z 

.  892 

.  107 

.956(-4) 

. 134(-3) 

Table  VI. 


Transaction  probabilities  |0,0>  -►  ( n ,  m )  j  j  '  ,  m  '  >  ,  and 
incoherent  and  diffractive  summed  results  for  an 
incident  beam  of  . 1  eV  molecules,  using  parameter 
set  3,  with  e=l  and  ^  =.022. 


|0.0> 

12, 0> 

|2. 1> 

12, 2> 

(0,0) 

.766 

.  128 

,«-16 

~10 

.831(-5) 

(±1.0) 

.311(-1) 

.459(-2) 

.626(-5) 

, 130(-4) 

(0,±1) 

.137(-1) 

.267(-2) 

.  150(-4) 

.813(-4) 

(±1,±1) 

.449(-3) 

.778(-4) 

.563(-6) 

. 189(-5) 

Dm 

.799 

.  198 

.9(-4) 

.496(-3) 

z 

.  857 

.  143 

.448(-4) 

.204 (-3) 

